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Abstract. We discuss the possibility to launch an outflow from the close vicinity of a protoplanetary core considering a model 
scenario where the protoplanet surrounded by a circum-p/a/ietary accretion disk is located in a circum-.s/e//ar disk. For the 
circum-planetary disk accretion rate we assume A/ cp ~ 6 x 10 Afjyr - implying peak disk temperatures of about 2000 K. 
The estimated disk ionization degree and Reynolds number allow for a sufficient coupling between the disk matter and the 
magnetic field. We find that the surface magnetic field strength of the protoplanet is probably not more than 10 G, indicating 
that the global planetary magnetosphere is dominated by the circum-planetary disk magnetic field of ^50 G. The existence 
of a gap between circum-planetary disk and planet seems to be unlikely. The estimated field strength and mass flow rates 
allow for asymptotic outflow velocities of ^60kms _ . The overall outflow geometry will be governed by the orbital radius, 
resembling a hollow tube or cone perpendicular the disk. The length of the outflow built up during one orbital period is about 
100 AU, depending on the outflow velocity. Outflows from circum-planetary disks may be visible in shock excited emission 
lines along a tube of diameter of the orbital radius and thickness of about 100 protoplanetary radii. We derive particle densities 
of 3000 cm -3 in this layer. Energetically, protoplanetary outflows cannot survive the interaction with a protostellar outflow. 
Due to the efficient angular momentum removal by the outflow, we expect the protoplanetary outflow to influence the early 
planet angular momentum evolution. If this is true, planets which have produced an outflow in earlier times will rotate slower at 
later times. The mass evolution of the planet is, however, hardly affected as the outflow mass loss rate will be small compared 
to the mass accumulated by the protoplanetary core. 
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systems: protoplanetary disks 



1. Introduction 

With the discovery of extrasolar planets during the last decade 
the scientific interest in planet formation has been increased 
substantially. In particular, with the help of the computational 
power existing today, it has become possible to perform nu- 
merical simulations of a circumstellar accretion disk contain- 
ing and building up an orbiting protoplanetary core (e.g. Kley 
[T999l Bryden e t al.[l999l Lubow et al. [l999l Kley e t al. l200Tl 
D'Angelo et al. l2"0"0"2l Tanigawa & Watanabe l^tMl D'Angelo 
et al. l2003all20Q3bl Bate et al. l2003l . Although differing in cer- 
tain aspects (spatial resolution of the region close to the proto- 
planet, number of dimensions treated) these simulations have 
provided us with the same general results. In particular, the 
simulations show how tidal interaction between the protoplanet 
and the disk material opens up a gap in the circum-stellar accre- 
tion disk along the orbit of the planetary core. Mass accretion 
in the circum-stellar disk, however, continues across the gap. 
The circum-stellar disk material entering the Roche lobe of the 
protoplanet becomes captured and finally accreted by the pro- 
toplanetary core. The simulations also demonstrate how the ac- 
cretion stream initially affected by strong shocks waves eventu- 



ally builds up a circum-planetary sub-disk in almost Keplerian 
rotation close to the planet. 

On the other hand, in the context of astrophysical jets 
it is well established that outflow formation is causally con- 
nected to the presence of an accretion disk and strong mag- 
netic fields (see Blandford & Payne 1982 Pudritz & Norman 
[19861 Camenzind [T990l Shu et al. fT994l . This statement holds 
for a wide range in outflow energy and spatial scale - for 
young stellar objects, microquasars and active galactic nuclei. 
Observations as well as theoretical models strongly suggest 
that astrophysical sources of outflows are (highly) magnetized. 
The magnetic field is responsible for acceleration and collima- 
tion and for lifting the matter from the disk into the outflow. 

In the scenario of planet formation within a circum-stellar 
accretion disk, the presence of a magnetic field can be expected 
as well. A central protostellar dipolar magnetic field of 1000G 
surface field strength may provide only about 20/iG at 5 AU 
distance. This (weak) field, however, may act as a seed field 
for dynamo action taking place in the circum-stellar accretion 
disk, in the protoplanetary core, or in the circum-planetary disk. 
Most of the present-day Solar system planets carry a substantial 
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magnetic field. We expect these fields to be present also during 
the phase of planet formation. 

Considering (i) the numerically established existence of 
accretion disks around protoplanetary cores as a natural fea- 
ture during the planet formation phase and (ii) the feasibility 
of a large-scale magnetic field in the protoplanetary environ- 
ment, the question arises whether the combination of these two 
factors may lead to the launch of an outflow similar to the 
phenomenon observed on much larger spatial and energetic 
scales. Up to date, the literature on this topic is rather rare. 
To our knowledge, the possible outflow activity from Jupiter- 
sized protoplanets has been investigated only by a single pa- 
per (Quillen & Trilling [L993 . Applying stationary models of 
magnetohydrodynamic outflow formation developed for prote- 
stors, these authors came to the conclusion that protoplanetary 
outflows might indeed exist with velocities ~ 20kms~ 1 and 
mass loss rates ~ 10~ 8 Jupiter masses per year. 

In the present paper we further investigate the feasibility 
of launching outflows from accretion disks around protoplan- 
ets. The paper is organized as follows. Section 2 summarizes 
the basic ideas of magnetohydrodynamic ally driven outflows 
and their possible application to protoplanets. We then discuss 
the circum-planetary accretion disk. Section 3.1 is devoted to 
numerical simulations of circum-planetary disk in the litera- 
ture. We estimate the temperature (Sect. 3.2) and ionization 
state (Sect. 3.3) of such disks. The next section considers the 
magnetic field of protoplanet (Sect. 4.1), and circum-planetary 
disk (Sect. 4.2), and the expected MHD properties of the out- 
flow (Sect. 4.4, 4.5). Based on the estimated properties of the 
system we discuss the parameter space for MHD driven out- 
flows from circum-planetary disks and their observational ap- 
pearance (Sect. 5). We conclude the paper with our summary 
(Sect. 6). 

We will consider a complex model scenario consisting of a 
number of constituents. In order to label their parameters, we 
use the following notation. Parameters of the planetary core 
have the subscript "p", those of the outflow the subscript "out", 
those of the circum-planetary disk the subscript "cp", those of 
the central star the subscript Parameters of the present- 
day Jupiter are denoted by the subscript "J". Further subscripts 
will be defined in the text. Variables without subscript denote 
general quantities. 

2. Magnetohydrodynamic (MHD) formation of 
outflows 

2. 1. The MHD model of astrophysical outflow 
formation 

We first outline briefly the essential features of MHD formation 
of outflows. The generally accepted "standard" model under- 
stands the outflow as magnetically driven (Blandford & Payne 
[T982I Pudritz & Norman fl986l Camenzind fT990l Shu et al. 
I1994I . Following current theoretical models the outflow origi- 
nates in the innermost part of a magnetized "star"-disk system 
(the "star" can be a young stellar object or a collapsed object). 
Whether the outflow magnetic field is primarily anchored in 
the accretion disk or in the surface of the central object, is not 



yet clear. The fact that outflows are observed also from sources 
containing a central black hole implies that the key magnetic 
field component is indeed provided by the surrounding accre- 
tion disk. The existence of non-relativistic outflows implies that 
the relativistic character of the source of an outflow cannot be 
a major constraint for its launching. 

We summarize the MHD formation of outflows process as 
follows. 

- The magnetic field necessary for launching can be gener- 
ated by dynamo action in the "star"-disk system or advected 
from the interstellar medium. 

- A fraction of the accreting plasma is lifted up to higher al- 
titudes by magnetic forces. The matter couples to the large- 
scale magnetic field, launching a slow wind. 

- The matter is then accelerated, first by magneto-centrifugal 
forces, then, further out, by Lorentz forces. Poynting flux is 
converted to kinetic energy. 

- Inertial forces lead to a "bending" of the poloidal field lines 
(i.e. the induction of a toroidal magnetic field). The toroidal 
field tension collimates the outflow into a narrow jet. 

- The plasma velocity subsequently exceeds the speed of the 
MHD waves. The asymptotic, super-magnetosonic flow is 
ballistic and causally decoupled from outer boundary con- 
ditions. 

- At the point where the outflow front meets the interstellar 
medium, a bow shock develops. 

It has been discussed in the literature that outflows can only 
be formed in a geometrical configuration underlying a cer- 
tain degree of axisymmetry (Fendt & Zinnecker ll998> . In the 
case of a circumplanetary disk-outflow-system, an axisymmet- 
ric alignment is indicated as the strong differential rotation of 
the circum-planetary disk would control the dynamo mecha- 
nism to provide magnetic fields (in the disk or in the planet) 
aligned with the rotational axis. 

Recent numerical simulations confirm this scenario of out- 
flow formation which has been developed over many years 
mainly by stationary MHD models. Axisymmetric MHD simu- 
lations have shown the self-collimating property of MHD flows 
from rotating disks (Ouyed & Pudritz 1997 ) also under the de- 
collimating effect of a turbulent magnetic diffusivity (Fendt & 
Cemeljic 2002 1. A stellar wind dominated disk-outflow does 
not collimate on the spatial scales considered in the simula- 
tions (Fendt & Elstner 20001 . Three-dimensional simulations 
of collimating disk winds prove that strong non-axisymmetric 
perturbations at the base of the outflow indeed lead to grow- 
ing MHD instabilities eventually disrupting the outflow flow 
(Ouyed et al. l2003l . 

2.2. The question of protoplanetary outflows 

In this section we identify the critical points for MHD out- 
flow formation in the context of protoplanetary outflows. First, 
one has to discuss the existence of a circum-planetary accre- 
tion disk and their parameters (accretion rate, temperature and 
degree of ionization). Circum-planetary accretion disks have 
not yet been investigated in detail theoretically or numerically. 
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Therefore, in our estimates we will apply the following ap- 
proach. We expect the outflow to be launched in the inner part 
of the circum-planetary disk. For this inner sub-disk, we as- 
sume a mass flow rate as determined by the numerical simu- 
lations of the circum- stellar accretion disk, and treat the disk 
structure assuming a standard, thin Shakura-Sunyaev disk. 

In the second step we will derived some estimates on the 
magnetic field structure close to the protoplanetary core. Two 
scenarios are plausible - an accretion disk equipartition mag- 
netic field or the field of the protoplanet itself. 

In general, the numerical simulations of the accretion 
stream around the protoplanet do not show any indication for 
outflow activity around the protoplanet. However, we know that 
astrophysical outflow are launched magnetically. Therefore, we 
cannot expect at all to detect such a process in a purely hydro- 
dynamic treatment, which, in addition, does not resolve the re- 
gion close to the planet. Radiation driven or gas pressure winds 
can be launched from a very hot accretion disk, however, such a 
scenario is unlikely for protoplanets where relatively cool disk 
temperatures are required to allow for solid dust grain conden- 
sation (Bryden et al. H999i . 

3. The circum-planetary accretion disk 

3. 1. Numerical evidence for circum-planetary 
accretion disks 

During the last years, a number of papers have been pub- 
lished considering time-dependent hydrodynamic simulations 
of a circum-stellar accretion disk which is enclosing an orbit- 
ing protoplanetary core (Bryden et al. 1999 , Lubow et al. 1999 
Nelson et al. l2000l Kley et al. l200"Tl D' Angelo et al. l2002l Bate 
et al. 120031 . In summary, all these simulations have demon- 
strated how the tidal interaction between the protoplanet and 
the viscous circum-stellar accretion disk opens a gap along 
the orbit. The feasibility of such a scenario has been predicted 
some decades ago (Lin & Papaloizou 1986 1993 ; Artymowicz 
& Lubow 1996 1. A similar dynamical evolution is known from 
simulations of circum-binary accretion disks (Artymowicz & 
Lubow [T994l 

Essentially, accretion of matter within the circum-stellar 
disk continues across the gap. That part of the material which 
spirals around the Roche lobe of the planetary core will then 
enter the Roche lobe and finally be captured by the planetary 
gravitational potential. The spiraling matter generates a density 
wave propagating into the circum-stellar disk. The path of the 
material approaching the Roche lobe follows more and more 
a circum-planetary orbital motion representing the feature of 
a circum-planetary accretion disk. Depending on the mass of 
the planetary core the gas orbiting within the circum-planetary 
disk reaches different rotational regimes. The simulations show 
that the orbital velocity profile can deviate from a Keplerian 
rotation by up to 10% primarily due to the smoothing of the 
gravitational potential (D' Angelo et al 2002 1. An almost circu- 
lar orbital symmetry is reached close to the central planet. On 
the other hand, at smaller radii the deviation from the Keplerian 
velocity is larger. 
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The governing parameter for our studies is the accre- 
tion rate of the circum-planetary disk. In general, we find a 
good agreement in the above mentioned papers concerning this 
value. Kley et al. ( 120011 give an accretion rate towards the plan- 
etary core of 

M cp ~ 6 x lO-Sfcfjyr- 1 , (1) 

a value which is confirmed also by three-dimensional simula- 
tions (D' Angelo et al. 2003a). The exact value depends on cer- 
tain model assumptions for the circum-stellar disk as the disk 
mass distributed along the orbit (3.5 x 10 _3 Mq between 2 
and 13 AU), the disk height (h(r) = 0.05 r) and the viscosity 
parameter (a~4x 10~ 3 ). Lubow et al. ©99 1 give accretion 
rates typical for their simulations of A/ cp ~ 4.5 x 10 _5 Mjyr _1 
for similar disk parameters as Kley et al. (1200 1> . Bryden et 
al. (11999 ) obtain accretion rates between M cp ~ 10 -5 Mjyr -1 
and M cp ~ 3 x 10 _4 Afjyr _1 which are, as expected, larger in 
the case of a higher disk viscosity. In the present paper, we will 
in general refer to the model parameters of Kley et al. (1200 1> 
and D' Angelo et al. (2002 ) and apply their numerically derived 
accretion rate of 6 x 10 -5 Mjyr -1 . This value is substantially 
higher than the value estimated by Quillen & Trilling dl998> . 

A difficulty with applying these numerical simulations 
to the circum-planetary accretion flow is their lack of suf- 
ficient numerical resolution in this region. In particular, this 
holds for the innermost part of the circum-planetary disk 
where we expect the outflow to be launched. The interest- 
ing length scales here are the planetary gravitational radius 
r H = L> p ( m p/ 3M *) 1/3 (Hill radius), 

the radius of the planetary core of presumably several Rj, and 
the orbital radius D p . 

Recent simulations by Tanigawa & Watanabe (2002), par- 
ticularly considering the gas accretion flow close to the pro- 
toplanet, apply a numerical grid extending from 0.025 rn to 
12 th resolving 0.005 rn close to the planet. For parameters 
similar to Eq. [2] this corresponds to about ~ 20i? j and cov- 
ers reasonably well the region where we expect outflow for- 
mation to happen. Due to a different setup - no gap in the 
circum-stellar disk along the orbit - the accretion rate towards 
the planetary core is higher by factor of hundred compared 
to the literature (Lubow et al. 119991 D' Angelo et al. 120021 . 
Compensating for this effect the authors confirm values for the 
re-scaled surface densities and accretion rates similar to Eq.0 
Further, Tanigawa & Watanabe deal with an isothermal disk, 
and, thus, cannot deliver the temperature profile in the disk. 
The authors claim that the main features of the flow dynamics 
are roughly the same for adiabatic and isothermal simulations 
(with a factor 3 difference in the accretion rate). 

Due to lack of numerical resolution close to the planet, the 
simulations published so far give little information about the 
structure of the circum-planetary disk (e.g. Kley 2002- Three- 
dimensional simulations are essential to determine the scale 
height (h/r) of the sub-disk. Further, the disk structure also de- 
pends on the number value for the viscosity parameter which is 
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generally taken constant throughout the whole computational 
domain (in particular, it is assumed to be the same for circum- 
stellar and circum-planetary disk). So far, the simulations have 
shown that the scale height of the sub-disk is much smaller 
compared to the main disk due to the gravitational potential of 
the planet (Bate et al. l2003l . 

Preliminary results of very high resolution numerical sim- 
ulations (Ciecielag et al. 2000 1 suggest that the spiral density 
wave launched in the outer regions of the circum-planetary disk 
may propagate into the inner disk regions. However, simula- 
tions by Bate et al. 12003) indicate that the strong spiral shocks 
observed in 2D simulations are greatly diminished in 3D simu- 
lations. In recent simulations by D'Angelo (2003 J the circum- 
planetary disk becomes almost axisymmetric in the region 
within 20 i?j from the planet. In the end, it is up to future high- 
resolution circum-planetary disk simulations to show whether 
the simple assumption of an axisymmetric, thin Keplerian disk 
must be revised. 

In summary, numerical hydrodynamic simulations of 
circum- stellar accretion disks give clear evidence for the ex- 
istence of circum-planetary disks in almost Keplerian rotation. 
In order to understand further the accretion process towards the 
protoplanetary core, the numerical simulation of the circum- 
planetary accretion disk itself is essential. No such work has 
yet been presented in the literature. Therefore, we have to rely 
on estimates obtained by scaling the standard accretion disk pa- 
rameters to protoplanetary parameters taking into account the 
results of the above mentioned circum-stellar disk simulations 
as model constraints. In the following we will consider a thin 
disk with accretion rate as given in Eq.^ 

3.2. Circum-planetary disk temperature and density 

Only matter of a sufficiently hot and, thus, sufficiently ion- 
ized accretion disk can couple to and interact with a magnetic 
field, eventually leading to the launch of winds or outflows. 
Therefore, we first need an estimate for the disk temperature 
profile. As discussed above, in order to estimate the tempera- 
ture profile of the circum-planetary accretion disk, we will con- 
sider the limit of a standard (thin, Shakura-Sunyaev) disk (see 
also Nelson & Benz 2003 1 with an accretion rate Eq.Hprovided 
by the numerical circum-stellar disk simulations. 

A first estimate of the circum-planetary disk temperature 
can be derived from the well known equation for the effective 
surface temperature of a geometrically thin disk surrounding a 
central body of mass M and an accretion rate M as a function 
of radius r is 

/ • \ 1/4 

3 GMM \ 

with the Stefan-Boltzmann constant a (see e.g. Pringle [T98U . 
For Jupiter-like parameters and a circum-planetary disk accre- 
tion rate as in Eq.^we obtain 

r, W = 1280K ( 6xW M s Z iyl J 0) 



Interestingly, we see that the "power" of the disk ~ MM 
normalized to the mass of the central object is similar for 
the circum-planetary and the circum-stellar case. However, the 
circum-planetary disk extends a factor (Rq/Rj) closer to the 
origin of the gravitational potential. Therefore, a high temper- 
ature and luminosity can be expected in the inner disk. 

The disk temperature as estimated in Eq.[3]is substantially 
larger than that of Quillen & Trilling d!998i . This is due to 
the much higher accretion rate we consider as motivated by the 
numerical simulations. The temperature derived in Eq. [3] is a 
measure for the disk surface temperature. 

Viscous heating of the disk will lead to even higher tem- 
peratures. For comparison, we show the mid-plane temperature 
distribution in the case of the a-parametrization for the viscos- 
ity (Shakura & Sunvaev ll973t in the appendix. As a result, we 
find a remarkable similarity between the circum-planetary and 
the circum-stellar accretion disks. 

A straightforward estimate of the mid-plane disk tempera- 
ture can be made by assuming radiative diffusion in the disk. 
Then the effective black body surface temperature is related 
to the central mid-plane temperature by = (3/8)rT c 4 ff 
(Hubenv I1990i . This gives mid-plane temperatures ~ 2 — 4 
higher than the surface temperature and rather insensitive to 
the vertical optical depth r. Our simple estimates yielding a 
rather high circum-planetary accretion disk temperature is in 
agreement also with disk models discussing the formation of 
the Galilean satellites (Canup & Ward l2002> . 

In order to derive a realistic disk internal temperature 
profile, more sophisticated theoretical model calculations are 
needed taking into account radiation transfer and a realistic de- 
scription of the opacity. Only few such models are available for 
protostellar accretion disks (e.g. D'Alessio et al. 119981 Malbet 
et al. l200T1 Dullemond et al. l2002l . none for circum-planetary 
disks. For the example of TTauri star accretion disks, it has 
been found (Malbet et al. 1200 1> that the central disk tempera- 
ture is a factor ten higher than the surface temperature within a 
disk radius of about 10 AU. Malbet et al. have pointed out that 
the other important heating source for the disk besides viscos- 
ity, namely the radiation of the central 4000 K star has not been 
taking into account in their model. D'Alessio et al. ( 119981 con- 
sider stellar irradiation as the main heating source of the disk 
upper layers. Similarly, Tanigawa & Watanabe (2002J have ar- 
gued that the temperature profile of a circum-planetary sub- 
disk within the Hill radius is substantially affected by heating 
of the protoplanetary core, in particular in the late stages of 
planetary growth. For our estimates, we do not consider the de- 
tailed vertical structure of the disk and the effect of irradiation. 
Our simple temperature estimate must be therefore considered 
as a lower limit. 

Recent 2D thermo-hydrodynamic simulations (D'Angelo 
et al. 2003b I were able to resolve the temperature profile close 
to the planet with reasonable resolution. These simulations take 
into account radiative heating and cooling as well as a sophis- 
ticated opacity model. The authors find a peak mid-plane tem- 
perature of 1200 - 1600K at radii of r = 5 - 10i?j from the 
protoplanet (D'Angelo 2003). Essentially, these temperatures 
are somewhat below the values given by our simple ansatz of a 
standard disk. This agrees with models of circum-stellar disks 
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taking into account radiative transfer and which give in gen- 
eral a lower mid-plane temperatures compared to the treatment 
without radiative transfer (Dullemond et al. I2002l >. 

A very high disk temperature might be problematic as po- 
tentially leading to a disk evaporation as discussed for disks 
around collapsed objects (e.g. Meyer & Meyer-Hofmeister 
119941 . Disk evaporation would affect just that (innermost) part 
of the circum-planetary disk where we expect the origin of the 
outflow. However, due to the fact that the circum-planetary disk 
reaches deeper into the gravitational potential of the central 
body, its a-disk scale height is less compared to the case of 
a circum-stellar disk (see eq. IA.5I >. This, and also the general 
similarity to the circum-stellar disk, indicates that the circum- 
planetary disk does not blow up and remains a thin disk close 
to the planet. 

3.3. The circum-planetary disk ionization 

A minimum degree of ionization is necessary to couple the 
matter to the magnetic field sufficiently. In the case of proto- 
stellar accretion disks several sources of disk ionization have 
been identified, (i) thermal ionization, (ii) ionization by cosmic 
ray particles, and (iii) ionization by non-thermal X-rays. As 
discussed above, the temperature and density profiles are com- 
parable in circum-stellar and circum-planetary accretion disks. 
Since the circum-planetary disk is located within the circum- 
stellar disk, it is straight forward to assume that also the ioniza- 
tion sources are the same, indicating that also the the degree of 
disk ionization will be similar in both cases. 

To determine the ionization degree of a protostellar accre- 
tion disk is a complex enterprise and depends on a number 
of unknown quantities (see e.g. Gammie 119961 D'Alessio et 
al. fT998l Fromang et al. l2002l l. In our case, we are interested 
in the innermost part of the circum-planetary accretion disk, 
where the launch of an outflow may take place. In this region 
the disk is hot and dominated by thermal ionization. The degree 
of ionization is then 

„ / a \V2 / T c \ 3/4 

t n n \ -i/2 exp(-25188T- 1 ) 
' U0 15 cm- 3 y 1.15 x 10-u ( ' 

with the neutral particle density n n and a the abundance of K 
atoms relative to Hydrogen (see Fromang et al.[2002|and ref- 
erences therein). Thermal ionization is dominant for mid-plane 
temperatures >10 3 K. Applying our previous estimates for the 
temperature profile in circum-planetary disks, such tempera- 
tures can be expected for radii r 65 Rj. Considering density 
values as given in Eq. (A. 2) we find an ionization degree above 
lCT 11 forr <65i?j. 

Note that the degree of thermal ionization of the circum- 
stellar disk at the orbital radius of 5.2 AU is much lower (Eq.|4ji 
and negligible. However, cosmic ray ionization may still pro- 
vide ionization levels up to 10~ n at these radii (Gammie 1 

1 There are still large uncertainties in the ionization rate along the 
disk surface, mostly connected to the possible impact of the stellar 
wind 



1996 see Eq. |5} and may therefore also contribute to the 
circum-planetary disk ionization. The vertical profile of the 
ionization level can only be calculated numerically. Model sim- 
ulations of circum-stellar disks indicate that lower disk layers 
may actually be decoupled from the magnetic field (Fromang et 
al. l2002> . This has already been discussed by Gammie ( 1996 > 
who pointed out sufficient coupling between matter and field is 
essential in order to excite the disk turbulence via the magneto- 
rotational instability (Balbus & Hawlev ll991> . The results by 
Gammie further indicate a two-layered accretion structure in 
protostellar accretion disks for radii > 0.1 AU. Regions in 
which the magnetic Reynolds number is ReM < 1 are mag- 
netically "inactive". 

It is interesting to work out these numbers for a circum- 
planetary disk. The magnetic Reynolds number is related to a 
certain degree of ionization. For the case of a circum-planetary 
disk, we obtain 

Re u = 7.6 x 10 15 £a 1/2 ( — V'Y -^—) ( —) ^ (5) 

Using the number values derived above, we see that in the inner 
part of the disk a sufficient ionization level £ <;10 -15 is always 
guaranteed providing a good coupling of matter and magnetic 
field as an essential condition in order to launch an outflow. 

4. The magnetic environment of the protoplanet 

We know from other astrophysical sources that outflow for- 
mation is basically a magnetic phenomenon - outflows are 
launched, accelerated and collimated by magnetic forces. 
Similar to the case of protostars and protostellar outflows, two 
model scenarios for the magnetic field structure around the pro- 
toplanetary core seem to be feasible. One is that of a magnetic 
field distribution dominated by the protoplanetary magneto- 
sphere. In the other, the accretion disk magnetic field is the 
main driver for the outflow. In the following we will elabo- 
rate in greater detail the scenario of a magnetohydrodynamic 
interrelation between the protoplanetary core, surrounding ac- 
cretion disk and a possible outflow. 

4. 1. The protoplanetary magnetosphere 

Many of the solar system planets carry a substantial magnetic 
field, in particular the outer gas giants. Jupiter has a well or- 
dered, strong dipolar magnetic field of about 10 G aligned with 
the rotational axis. It must be mentioned, however, that the de- 
tailed structure and orientation of the solar system planets dif- 
fers substantially and the theoretical understanding of the un- 
derlying (necessarily diverse) dynamo processes is not yet un- 
derstood (e.g. Connernev l 19931 Moss & Brandenburg l 1 995V 

So far, we are not aware of any theoretical treatment about 
the magnetic field evolution in protoplanets. Therefore, our 
considerations of the structure and field strength of the proto- 
planetary magnetosphere have to rely on simplifying estimates. 

Nevertheless, we know from the comparison of stars and 
protostars that the large scale magnetic field strength observed 
in young stars is generally higher than the field strength found 
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in main sequence stars. The main reason is that the strength 
of stellar dynamo action depends on convection (turbulence) 
and rotation (see Bouvier 119901 Ghosh 119951 and references 
therein). These parameters are clearly more pronounced in the 
early phase of star formation. It might be straight forward to 
assume a similar extrapolation for the early formation phase 
of a planet. This would imply that the large-scale protoplane- 
tary magnetic field has been larger than that of a present-day 
planet. The seed magnetic field for the protoplanetary dynamo 
can be provided by the magnetosphere of the central star and/or 
the circum-stellar accretion disk where the planet is embedded 
in. What kind of dynamo (a 2 or ail) is operating in a proto- 
planet is not known. Similar to protostars a dynamo driven by 
turbulent convection (a) is therefore feasible. Whether, in ad- 
dition, differential rotation (Q) plays a role remains unclear. It 
has been shown that fully convective protostars do not rotate 
differentially (see below). There is theoretical indication from 
the planet formation that a giant gaseous protoplanet reaches 
an almost fully convective state within some hundred years 
(Wuchterl et al. I2000i . However, as for protostars, the pres- 
ence of an accretion disk and the inferred magnetic coupling 
between disk and planet will supply a substantial amount of 
differential rotation. Thus, both kinds of dynamo action can be 
expected, an ail-dynamo producing an aligned axisymmetric 
dipole, or an a 2 -dynamo with a non-axisymmetric field struc- 
ture. Also a quadrupolar field structure are feasible, in particu- 
lar for high Taylor numbers (Grote et al. H999i . 

The maximum field strength which could be generated by 
a dynamo can be estimated considering the fact that convection 



must be super- Alfvenic (Camenzind 1 1 9971 see also Kiiker & 
Riidiger 1999 ). The Alfven speed for this critical (mean) field 
strength is va — B cr / 1 and has to be compared to the con- 
vective velocity defined by the convective turnover time scale 
T con . The convective velocity v c can be derived from the energy 
transport in the protoplanet. Considering the convective radia- 
tion flux in the limit of mixing length theory, the luminosity of 
the protoplanet is L p = e c AI p v^/Rp. 

Theoretical models of the planetary evolution have shown 
that the luminosity of a young planet can be very high, L p ~ 
10~ 4 Lq (Burrows et al. 1200 It . Therefore, we estimate the 
mean convective velocity in a protoplanet as 
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Comparison to the mean Alfven speed in the planet provides 
the critical magnetic field strength for a hypothetical protoplan- 
etary dynamo, 
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the field strength observed in present-day solar system plan- 
ets. If this similarity is true, it indicates a substantial difference 
between the magnetic evolution of protostars and protoplan- 
ets. Secondly, this field strength is well below the equipartition 
field strength in the circum-planetary accretion disk (see be- 
low). Therefore, we believe that, in difference to the case of 
protostars, the protoplanetary magnetic field cannot dominate 
the dynamics of the system even close to the planetary core. 

It has been found that a good indicator for the level of mag- 
netic activity of a star with rotational period P is the Rossby 



number Ro = P/r con (see Pizzolato et al. 20Q3JI. If we adapt 
this relation for the case of a protoplanet we find 



Ro = 2.3 x 10" 



Qk(Rf 



30 



1/3 



io-u 



M n 



-5/6 



Rr 



5/6 



1/3 



(8) 



(7) 



with the Alfven Mach number of the convective motion Mp = 
w con/^A- A higher field strength would imply a magnetic 
quenching mechanism suppressing the convective motion (and 
the turbulence) which drives the dynamo. 

This estimate of an upper limit for the magnetic field 
strength in fully convective protoplanets is interesting in two 
respects. Firstly, the derived field strength is not higher than 



lO" 3 M y \3Rj 

Compared to protostars (Camenzind ll997> . this value is some- 
what higher (factor 3/2), indicating that the protoplanet can ex- 
hibit a high magnetic activity as the central protostar. However, 
the high level of activity is probably compensated by the 
weaker energy output due to the relatively low surface mag- 
netic field strength. 

So far we have neglected the effect of dissipative processes 
for the dynamo mechanism. Diffusive and viscous effects will 
certainly determine the dynamical evolution of the system. 
These effects are widely unstudied, and probably not yet fully 
understood (see e.g. Blackman 120031 Vishniac et al. I2003i . 
Therefore, in the following we will discuss some general char- 
acteristics of the protostellar dynamo by estimating dimension- 
less dynamo numbers (see Camenzind ll997l Kiiker & Riidiger 
1999 1 in the framework of the mean field dynamo theory. 

The first point to note is that the low Rossby number 
Ro << 1 indicates the protoplanet as a fast rotator. Thus, the 
scale of the dynamo a-effect, the a-dynamo parameter ap>, is 
given by the convective velocity ap> — v c . Another interesting 
question is whether the a or il-effect plays the major role for 
the dynamo. This can be estimated by comparing the magnetic 
Reynolds numbers R a = ap,Rp/rj and Rq = AilpRp/r], 
where r\ is the turbulent magnetic diffusivity in the planet and 
Ai7p the differential rotation. As a matter of fact, the latter 
quantity is quite unknown. For AOp ~ Hp we find a maxi- 
mum ratio (Rn/Ra) — QpRp/v c ~ 10 2 . On the other hand, 
it has been shown that fully convective protostars rotate al- 
most rigidly (Kiiker & Rudiger ll997> . This holds in general for 
spherical stars without latitudinal gradients in pressure, density 
and temperature. We are attempted to extrapolate these results 
to fully convective protoplanetary cores. As a consequence, 
AOp < < Hp and the mean field dynamo is dominated by the 
aip/ia-effect (turbulence). 

The a-effect becomes more efficient for faster rotators, as 



noted by Kiiker & Riidiger dl999t for protostars. Also this 
statement can certainly be generalized from protostars to pro- 
toplanets. In difference to protostars, the protoplanetary core 
is colder. Thus, the magnetic Prandtl number (measuring vis- 
cosity in terms of diffusivity) is generally small. For planets 
we have only P m = {v/r}) ~ 10~ 5 ...10~ 2 (see also Kirk & 
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Stevenson 1 19871 Starchencko & Jones 2002 ). In turn, this im- 
plies that the time scale for the turbulent planetary dynamo is 
very short. We expect a rapid increase in field strength during 
the linear regime up to a substantial value, and later a slow sat- 
uration to the maximum field strength (non-linear quenching of 
the dynamo-driving turbulence). 

4.2. Magnetosphere-disk interaction 

We now address the question whether the planetary magnetic 
field may, similar to the case of protostars, disrupt the inner part 
of the disk forming a gap between the planet and an inner disk 
radius. The gas pressure p in a standard disk with an accretion 
rate M and a disk half thickness h(r) around a central body of 
mass M is (neglecting factors of orders of unity) 

, . 1 M [GM\ 1/2 (h\~ l 

(Camenzind 1 1 9901 see also Riidiger et al. [1995J. The quantity 
a is the usual disk viscosity parameter for the viscous shear 
stress t r <f, = —ap. For a dipolar magnetic field with B ~ r -3 , 
pressure equilibrium between the magnetic pressure of the pro- 
toplanetary magnetosphere and the accretion disk gas pressure 
defines the inner disk radius r ln . In units of the radius of the 
central object R, the inner disk radius is located at 

Here, B s is the magnetic field strength on the surface of the 
central body. If the magnetic pressure is strong enough, it will 
open up a gap between the surface of the central body and the 
disk. In the case of protostellar magnetospheres, the inner disk 
radius is located at about 3 protostellar radii (Camenzind 1990 
Konigi n"991l l. However, for the case of a circum-planetary ac- 
cretion disk and a 10 G central dipolar planetary magnetic field 
the situation would be quite different, 

n„ n11 2 / 7 f B S \ i/7 ( M cp Y 2/7 nu 

(My 1 ' 7 (Rp_\ 5/7 (h/r\ 2/7 

\Mj) \3RjJ \0.l) 

Such a low (academic) value for the inner disk radius indicates 
that for the accretion rate considered here the dynamical impact 
of the planetary magnetosphere for the disk structure is small. 
Only for a 50 G surface field strength in combination with an 
accretion rate 1 00 times lower the inner disk radius moves away 
from the surface of the planet to a radius r- in ~ 2Rp. This esti- 
mate is consistent with the results of Quillen & Trilling ( 1998 1 
who derive ~ 1.3i?p. Quillen & Trilling argue that during 
later evolutionary stages during the planet formation the accre- 
tion rate will decrease, allowing for a disk disruption. However, 
we like to point out that a lower accretion rate also implies a 
lower mass loss rate for the outflow. 

The previous estimates demonstrate again how critical the 
scenario of protoplanetary outflow formation depends on the 



circum-planetary accretion rate. More detailed numerical sim- 
ulations of the circum-planetary accretion disk evolution are 
essential in order to clarify this point. So far, we have to rely 
on the simulations presented in the literature which give a high 
accretion rate. 

Although we have shown that it is unlikely that a protoplan- 
etary magnetic field can open up a gap in the circum-planetary 
accretion disk it is interesting to calculate the co-rotation ra- 
dius r co of such a magnetosphere. For a rotational period of the 
protoplanet of P p = 10 h , we find 

( M \ 1/3 / P \ 2/3 
r co = 2.25i^J - (12) 

which is comparable to its radius. The ratio r co /Rj is about a 
factor 10 larger than for the protostellar case. Note that break- 
up rotational period for a Rp = 3 i?j (i?p = 2 Rj) protoplan- 
etary core is about 15 h (8 h). 

As shown above, the circum-planetary disk most probably 
reaches below the co-rotation radius, implying an angular ve- 
locity at the basis of an outflow higher than the one of the cen- 
tral planet - in difference to the case of protostars where the 
outflow emanates from regions outside the co-rotation radius. 

It has been claimed in the literature that a stellar magnetic 
field sufficiently strong to open up a gap between the star and 
the disk would be automatically strong enough (e.g. Quillen & 
Trilling ) to launch a magnetized wind. However, one has to 
accept the fact that there is (yet) no numerical proof of such 
a statement. Stationary state solutions cannot trace the evolu- 
tion of an outflow (Camenzind 119901 Shu et al. 119941 Fendt 
et al. I1995i and time-dependent simulations of outflow forma- 
tion either do not consider a dipolar central field (e.g. Ouyed 
& Pudritz 119971 or do not show evidence for a collimated 
outflo w (Mill er & Stone [T9971 Goodson et al. [T997l Fendt & 
Elstner l2000> . The fact that strongly magnetized neutron stars 
and white dwarfs do not have collimated outflows is another in- 
teresting feature in this respect. On the other hand, outflows are 
observed in active galactic nuclei emanating from the accretion 
disk around a black hole which does not generate a magnetic 
field at all. 

4.3. The circum-planetary disk magnetic field 

Circum-stellar accretion disks are believed to carry a substan- 
tial magnetic field. Under certain circumstances, this disk field 
can act as a driver for the observed outflow activity. Further, a 
disk magnetic field of moderate strength plays the essential role 
of generating the disk turbulence by the magneto-rotational in- 
stability (Balbus & Hawlev Tl991> - which enables the disk 
angular momentum transport and eventually allows for accre- 
tion of matter at all. In turn, turbulent motion in a rotating ion- 
ized medium may give rise to "mean field" dynamo action (a- 
effect). 

As discussed above (see Sect. l3.2l and l3.3t there exist strong 
similarities between circum-stellar and circum-planetary accre- 
tion disks concerning the degree of ionization and the tempera- 
ture distribution. In the following we will assume that dynamo 
action in a circum-planetary disk is indeed feasible. Another 
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possibility to accumulate magnetic flux in the circum-planetary 
disk is by advection of ambient field (i.e. the magnetic field in- 
duced in the surrounding circum-stellar disk) within the accre- 
tion stream. 

From Eq.|9]we derive a circum-planetary disk equipartition 
field strength of 



B eq (r) = 480Ga~ 1/2 
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This field strength is definitely higher than the surface field 
strength estimated above for a protoplanetary dynamo, in par- 
ticular, since we have to consider a turbulence parameter a < 
1. We may therefore conclude that a hypothetical protoplane- 
tary outflow is launched from the disk and not from the planet 
itself. The actual disk magnetic field strength is certainly be- 
low the equipartition field and depends on the detailed proper- 
ties of the dynamo (see e.g. Riidiger et al. I1995t . In general, 
the toroidal disk magnetic field strength below the equiparti- 
tion value and the poloidal field strength one or two orders of 
magnitude below the toroidal field, with B r ~ B z depending 
on the dynamo number. Adding up, this gives a estimate for 
the circum-planetary accretion disk poloidal magnetic field at 
r = 10 Rj of about 10 - 50 G, if we consider a ~ 0.01. 

A disk field strength similar to the equipartition field 
strength would lead to disk instabilities. The disk magnetic 
field amplified by the magneto-rotational instability induced 
turbulence will be saturated by the buoyant escape in vertical 
direction leading to the formation of a magnetized disk corona 



(Merloni 2003). Further, disks with high Poynting flux imply- 
ing a relatively strong poloidal magnetic field are favored by 
a low viscosity. These results obtained for disks around stellar 
mass black holes may be tentatively generalized to other pa- 
rameters regimes. Three-dimensional MHD simulations show 
that a magnetic torus with initially equipartition field strength 
will change into a state in which the global field strength is 
(1/5) of the equipartition field strength (Machida et al. I2000i . 
However, locally, magnetic pressure dominated filaments can 
be found. 

If the differential rotation of the disk (fi-effect) completely 
dominates the magnetic field induction by the dynamo, the 
magnetic disk flux will be negligible, and the outflow might be 
launched as a magnetic "plasma gun" (see Contopoulos l 1 9951 . 
If the a-effect generates a substantial poloidal field component, 
we can expect an outflow to be launched similar to the well- 
known magneto-centrifugal mechanism (Blandford & Payne 
[19821 . 

4.4. Magnetization and asymptotic outflow velocity 

The essential parameter determining the asymptotic speed of 
the MHD outflow, is its magnetization, 



$ 2 ft| 



(Michel 1969|. Here, fip is the angular velocity at the foot 
point radius tf of a magnetic field line, $ = r F Bp refers 
to the magnetic flux from that area and M out to the outflow 
mass loss rate. For an outflow expanding in spherically ra- 
dial direction, Michel (1959) derived an analytical relation be- 
tween the asymptotic flow velocity follows and the magneti- 
zation, Vac ~ cr 1 / 3 (Michel scaling). We emphasize the im- 
portance of having a rapid rotation (51f) in order to launch a 
fast outflow. Faster rotators launch faster outflows. It has been 
shown that the power law interrelation above is maintained 
also forcollimating outflows, however with a different power 
index. In fact, the Michel scaling can be applied for a wide 
range of astrophysical outflows from protostellar outflows with 
v out ~ 500km/s and a ~ 10~ 8 to the extragalactic case with 
w out ~ c and a ~ 10 Fendt & Camenzind J1996i . 

Applying the Michel scaling for the asymptotic outflow ve- 
locity from circum-planetary accretion disks, we obtain 
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Here, we have considered the magnetic flux $ C p = 5 x 
10 22 Gcm 2 corresponding to a poloidal magnetic field com- 
ponent of B p ~ 10 G and a Keplerian period at r = 10i?j. 
For the outflow mass flow rate we assume 0. 1 % of the disk ac- 
cretion rate, a typical value derived in MHD outflow formation 
models. For comparison, the asymptotic speed of an outflow 
launched from the planetary surface would be ~ 13 km s _1 
only, with the magnetic flux $p = 5 x 10 20 G cm 2 and a rota- 
tional period of 10 h (present-day Jupiter). 

The asymptotic outflow velocity as derived in Eq. [factu- 
ally represents a lower limit. Taking into account a higher disk 
field strength and a more rapid disk rotation closer to the planet 
the asymptotic velocity can be larger by a substantial factor. 



4.5. Comparison to simulations of MHD outflow 
formation 

In general, two approaches have been undertaken to simulate 
the formation of an MHD outflow. One way is to treat the 
acceleration and collimation process above the disk and con- 
sider the disk "only" as a boundary condition for the mass 
flow and the magnetic field structure fixed in time. In the case 
of a monotonous disk magnetic field distribution a collimated 
outflow evolves out of a slow injection of matter in the disk 
corona (e.g. Ouyed & Pudritz IT9971 Fendt & Cemeljic |2UU2), 
The other approach includes the evolution of the disk in the 
simulation and the disk-outflow connection, i.e. the launching 
mechanism, can be investigated (Hayashi et al. 1996 Miller & 
Stone 119971 Kuwabara et al. I2000i . However, mainly due to 
numerical problems with the accretion disk stability, this ap- 
proach does not yet allow to cover the long-term evolution of 
the disk-outflow system. 

All these simulations are performed considering normal- 
ized variables. Having length scales normalized to an "inner 
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Fig. 1. Model scenario of hypothetical outflows from circum-planetary accretion disks. A protoplanet surrounded by a circum- 
planetary disk ("CP-disk") is located within the circum-stellar disk ("CS-disk"). Depending on the outflow velocity achieved, 
two different types of outflow activity is possible - a large scale outflow, possibly collimated, or a bipolar blob orbiting with the 
planet. 



disk radius", velocities to the Keplerian velocity at this ra- 
dius and densities corresponding to a certain, given mass flow 
rate (Ouyed & Pudritz I1997t . the resulting simulations may 
be scaled from protostellar outflow to the AGN, supposed that 
the derived parameters (as Mach number in the disk or outflow 
and disk temperature) are consistent with the observations. The 
question arises whether we can apply recent results of magne- 
tohydrodynamic simulations of outflow formation (e.g. Ouyed 
& Pudritz [T997l Fendt & Elstner l2000l Fendt & Cemeljic|2(f)2} 
also for the model scenario of a protoplanet. 

Essentially, there are three governing model assumptions 
involved in the simulations, which are quite well known from 
the observations. This is the fact that the accretion disk launch- 
ing the outflow is geometrically thin, implying that the sound 
speed is small compared to Keplerian speed. Further, the out- 
flow is "cold", i.e. gas pressure can be neglected (low plasma 
beta). The third assumption concerns the mass load of the out- 
flow. It is clear that we do not know these parameters for the 
case in the case of planetary outflows due to the (yet existing) 
lack of observational evidence. As a working hypothesis, we 
may assume so far that outflow formation works similarly on 
all scales and that the MHD simulations of protostellar jet for- 
mation can indeed be applied for the case of protoplanets. The 
main results we can extrapolate are then (i) asymptotic out- 
flow velocity several times the Keplerian velocity at the inner 
disk radius (Ouyed & Pudritz [T997i and (ii) a rough equipar- 
tition between magnetic and kinetic energy in the asymptotic 
flow. For protoplanetary outflows launched from a disk ra- 
dius of 10 Rj this would imply a asymptotic speed of at least 
15kms _1 . 



5. Observational appearance of protoplanetary 
outflows 

Finally, we discuss possible observational features of proto- 
planetary outflows based on the model scenario developed 
earlier in this paper. A Jovian mass protoplanetary core sur- 
rounded by a circum-planetary accretion disk orbits the central 
protostar along a Jupiter orbit. Tidal interaction between the 
protoplanet and the circum-stellar accretion disk has opened- 
up a gap with in the circum-stellar disk. A large scale magnetic 
field below the equipartition field strength is provided by the 
circum-planetary accretion disk. The outflow is launched mag- 
netically from the circum-planetary accretion disk. 

5. 1. Outflow kinematics and geometry 

A critical measure for the outflow velocity are the escape ve- 
locities. The escape velocity of a protoplanetary outflow from 
the gravitational potential of the protoplanet is 

/ r \ - V2 / M \ 1/2 
, e , j( r) = 59.6k mS -%J [ Wj ) ■ d6) 

This value is similar to the lower limit for the asymptotic flow 
velocity derived above for the Michel scaling, Eq.^] In com- 
parison, the escape speed for the protoplanetary outflow from 
the gravitational potential of the central star at the orbit of 
Jupiter is 

^(r) = 18.5kms- 1 (^) ^ ' (^f \ (17) 
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Therefore, a protoplanetary outflow propagating faster than the 
escape velocity from the protoplanet will also be able to leave 
the gravitational potential of the central star (see Fig. la). 

For a standard thin circum- stellar accretion disk of aspect 
ratio h/r ~ 0.1, the scale height at the orbit of Jupiter is h ~ 
0.5 AU ~ 210 Rj. During the early phases of gap formation, 
the width of that gap will set an upper limit for the radius of the 
protoplanetary outflow. If the outflow is launched during later 
stages of gap formation, MHD hoop-stresses may collimate the 
outflow to a radius smaller then the gap size. 

During one orbital period a protoplanetary outflow would 
propagate "length" of the outflow, is 



Zout.J = Voo-Porb,J = 150 AU 



60 kms 



(18) 



{"length" of the outflow). This is a substantial distance if we 
consider the size of the outflow origin of only several proto- 
planetary radii. If the life time of outflow is long enough (i.e. 
lasting more than several orbits), we expect the overall struc- 
ture of the outflow to be a hollow tube 2 perpendicular to the 
circum-stellar accretion disk with a diameter of the order of the 
orbital radius. The protoplanetary jet propagates within a layer 
along the tube/cone. The thickness of that layer is defined by 
the radius of the intrinsic outflow and is probably several hun- 
dred planetary radii, or about 1 % of the tube diameter. Here we 
see a major difference to protostellar jets where the mass flow 
is distributed over the whole width of the global structure. 

For an outflow velocity too low in order to leave the plane- 
tary gravitational potential, f CS c,j > v out, the material ejected 
from the disk will most probably be accumulated in a blob 
around the planet (see Fig. lb). As the blob remains bound to 
the planet, it will follow the orbit around the central star. The 
actual size of the blob will depend on the detailed balance be- 
tween ejection of material out of the circum-planetary disk and 
(spherical?) accretion of this matter back to the protoplanetary 
system. The maximum size of the blob is the Hill radius. The 
possible existence of such an extended sphere of hot, shocked 
gas around the protoplanet is quite interesting from the obser- 
vational point of view. 



The kinetic luminosity P^, 
tary outflow is about 
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(19) 



In comparison, for protostellar outflows with u out ~ 
500 kms - and M ut — 10 -8 Mq yr _1 this value is four or- 
ders of magnitude higher. We therefore conclude that a proto- 
planetary outflow could be detected only if launched during a 



Depending on the separation of the planet from the central star, 
the outflow geometry may be characterized better as a hollow cone. 
Centrifugal imbalance with the central protostar will open up the tube 
geometry into a cone as the outflow ascends from the ecliptic plane 
(see Fendt & Zinnecker 1998, Masciadri & Raga 2002 for a similar 
discussion in the context of outflows from binary protostars) 



time period when no protostellar outflow is present. Otherwise 
the stellar outflow launched in the inner part of the circum- 
stellar disk and expanding outwards to about 50 AU will just 
disrupt the narrow protoplanetary outflow. 

An essential property of magnetic winds and outflows is the 
ability to remove angular momentum very efficiently. This is 
due to the large lever arm, which is defined by the Alfen radius 
ta and typically a factor ten larger then the foot point radius of 
the corresponding magnetic field line. The angular momentum 
carried away per unit mass is 



L 



rBpB^ 



Q F r A 



9.5xl0 18 cm 2 s" 1 
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(20) 



where Slj? is the angular velocity at the foot point radius rp of 
the field line. The is a factor {ta/tf) higher then in the pure 
hydrodynamic case with L ~ Vlpr F . 

The total angular momentum flux J = ML carried away 
by the protoplanetary outflow can be estimated as 



3xl0 35 gcm 2 s~ 2 




With that, the total angular momentum loss over a orbital pe- 
riod at a Jupiter distance is about J out = 10 44 gcm 2 s~ 2 . 

In comparison, the (rotational) angular momentum of 
Jupiter is Jj = 6.3 x 10 45 g cm 2 s -2 (approximated by a homo- 
geneous sphere). Even if the mass carried away by the proto- 
planetary outflow is small compared to the mass of the central 
planet, the angular momentum of both components can be of 
similar order. This leads to the conclusion that a planetary out- 
flow launched for many orbital periods of the protoplanet may 
well affect the angular momentum evolution of the planetary 
core and, thus, the time scale for planet formation. 

5.2. Outflow density and observational features 

The time scale of the circum-planetary accretion disk can be 
roughly estimated to be of the order of 



r disk ~ Mj/M cp ~ 2 x 10 V ~ 1400P orb ,j. 



(22) 



Since the accretion rate decreases with time this value provides 
only a lower limit. It is clear, that the life time of the disk pro- 
vides a natural upper limit for the life time of the outflow 3 . 
The amount of mass deposited by the outflow into the ambient 
medium during its life time T ou t is 



M n 



M out r out ~ 10" 5 Mj 



10- 3 r d 




M n 



10" 2 M, 



(23) 
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" Statistical examinations as well as the kinematic evolution of pro- 
tostellar jets show that the life time of the outflow must be consider- 
ably shorter than the disk time scale 
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On the other hand, from the outflow mass accumulated dur- 
ing one orbital period being confined in a 100 i?j-wide sur- 
face layer along the global outflow tube of 10 AU diameter and 
150 AU length, we derive mean a particle density within this 
layer of about 2800 cm -3 . This value is surprisingly similar to 
protostellar outflows (see e.g. Mundt et al. I1990> . As the ve- 
locities are not too different from protostellar jets, we therefore 
expect the existence of shock excited optical forbidden lines 
also in protoplanetary outflow. The observed features would, 
however, look different. In protostellar outflows we observe a 
wide shock structure extending over the whole 50 AU width of 
the outflow. In comparison, in the case protoplanetary outflows 
one has to consider a bunch of "shocklets" distributed around 
the hollow tube of the flow. As the kinetic luminosity of the 
proposed protoplanetary outflows is at least two orders of mag- 
nitude lower than for protostellar outflows, we expect the same 
for luminosity in the shock emission. However, as discussed 
above, one can suspect that protostellar outflows are launched 
at times scales earlier than those considered the formation of 
protoplanets and their hypothetical outflows. 

5.3. Influence on the planet formation process 

Can protoplanetary outflows affect the process of planet forma- 
tion? Basically, the outflow restores accreted matter back into 
the ambient medium while efficiently removing angular mo- 
mentum from the accretion disk. 

Considering the number values derived above we come to 
the conclusion that the mass evolution of the planet cannot be 
affected by the mass loss of the outflow as it is only about 1% 
of the accretion rate. 

The impact of the enhanced angular momentum loss due to 
the magnetic outflow can be twofold. First, the time-scale for 
the accretion process can be decreased as the disk angular mo- 
mentum is removed more efficiently compared to purely vis- 
cous transport in the disk. The acceleration in the growth rate 
for the protoplanet is difficult to estimate. A first guess would 
be that this time-scale is lowered by the same factor of 100 
at which the angular momentum loss is increased. Secondly, 
the angular momentum loss by the disk wind will influence the 
planet rotational period. If angular momentum is redistributed 
into the interstellar medium and, thus, not accreted by the plan- 
etary core, the final planet can be expected to rotate slower. 

6. Summary 

In this paper we have considered the possibility of launching 
magnetized outflows from circum-planetary accretion disks. 
We discuss a model scenario where the protoplanet is accom- 
panied by a circum-planetary accretion disk which is fed from 
the surrounding circum-stellar accretion disk. This scenario is 
motivated by recent numerical simulations of planet formation. 

For the outflow formation itself we suppose a mechanism 
similar to other astrophysical outflows (protostellar jets, extra- 
galactic jets) where the flow is magnetically launched as a disk 
wind and then accelerated and collimated into a narrow beam 
by Lorentz forces. 



We have estimated the magnetic field strength of a fully 
convective protoplanetary core and find an upper limit for 
the surface field strength of about 10 G. For the surrounding 
circum-planetary accretion disk we assume a accretion rate 
of M cp = 6 x 10~ 5 Mj yr for Jupiter-mass protoplanet. The 
equipartition magnetic field strength in such a disk is of a few 
100 G and larger than the protoplanetary magnetic field. The 
magnetic flux from the circum-planetary accretion disk is esti- 
mated to 5 x 10 22 G cm 2 . 

We further investigated the circum-planetary accretion disk 
temperature, its ionization state and magnetic Reynolds num- 
ber. The circum-planetary disk temperature may reach values 
up to 2000 K. We find strong indication for a sufficient matter- 
field coupling underlining the magnetic character of the disk. 
The latter is an essential condition for outflow launching. 

From the estimated accretion rate and magnetic field 
strength we consider the outflow magnetization for the asymp- 
totic outflow velocity. Applying the Michel scaling for mag- 
netic outflows this velocity is about 63 km s _1 and of the order 
of the escape speed for the protoplanet. However, a modified 
Michel scaling implied by the collimated structure of the flow 
may result in even higher velocities. 

For reasonable estimates for the outflow mass loss rate of 
M ut = 10~ 2 M cp and a velocity of i> ou t = 60kms , the 
estimated kinetic power of a protoplanetary disk outflow is 
Pkin — 10 30 ergs _1 . These estimates also imply that during 
one orbit of the planet a mass of about 10 -5 Jupiter masses can 
be deposited in the interstellar space by the protoplanetary out- 
flow. In the same time period this outflow would travel about 
150 AU. 

Two outflow scenarios are feasible depending on the out- 
flow velocity. If the outflow cannot escape the protoplanetary 
gravitational potential, we expect the flow to build up a ex- 
tended (bipolar) blob of hot gas orbiting with the planet. If the 
outflow exceeds the escape speed of the protoplanet, a colli- 
mated outflow may be formed which (at a planetary orbit of 
about 5 AU) can also escape the stellar potential. 

The overall outflow geometry built up by the orbiting 
source of a fast outflow is that of a hollow tube or hollow cone 
perpendicular to both circum-planetary and circum-stellar disk. 
For the jet mass loss rate considered we find particle densities 
in the layer along the tube of the outflow of 2800 cm -3 , a value 
similar to the density in protostellar jets. Shock excited forbid- 
den line emission can therefore be expected. 

Energetically protoplanetary outflows cannot survive the 
interaction with a protostellar jet which is launched from the 
inner regions of the circum-stellar accretion disk and can, thus, 
be only present if the stellar outflow has ceased to exist. 

The efficient angular momentum loss by the magnetized 
disk wind may affect the accretion time scale and, thus, the 
time scale for planet formation. The angular momentum loss 
by the disk wind or outflow per one orbital period is about the 
rotational angular momentum of today's Jupiter. 

In summary, our model estimates strongly rely on the ac- 
cretion rate of the circum-planetary accretion disk. We ap- 
ply M cp = 6 x 10~ 5 Mj yr, a value which is supported 
by many independent numerical simulations of planet for- 
mation in accretion disks. Further model constraints can be 
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expected from future, high resolution numerical simulations 
of the circum-planetary accretion disk clarifying the detailed 
circum-planetary disk structure close to the protoplanet. 
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Appendix A: The standard disk parameters 



disk, 



Here we apply the Shakura & Sunyaev d!973> accretion disk 
parametrisation of a thin a-viscosity disk to the case of circum- 
planetary disks. We show the temperature and density profile 
and for comparison the same also for the circum-stellar disk. In 
our case radiation pressure can be neglected. Free-free absorp- 
tion dominates electron scattering. The sound speed is given by 
c s = {kBT/m pT0 t). The particle density is 



n(r) = 1.6 x 10 14 cm- 3 a- 7 / 10 



M 



11/20 



6 x lO-SAfjyr- 1 / 

M \ 5/8 / „ N -15/8 



Mi) 



2.8 x 10 14 cm- 3 a- 7 / 10 



15 Ri 



M 



(A.l) 



11/20 



1.2 x IO-TMqvt- 1 
M\ 5/8 / r ^ 15/8 



Mr. 



while the disk temperature is 



T(r) = 2800Ka~ 1/5 







M 



15 Ro 







3/10 



(A.2) 



6 x lO-SMjyr- 1 
M 

M 

1.2 x lO^MQyr- 1 



3700K a 



-1/5 



15i?,, / 

\ 3/10 



(A.3) 



A/\ 1/4 / r ^" 3/4 



Mr. 







15 Re 







(A.4) 



The temperature derived here is the mid-plane disk temperature 
and is in agreement with the simple estimate of the surface tem- 
perature in Eq. [3] but also the more sophisticated disk model 
probing the vertical circum-stellar disk structure (D' Alessio et 
al. H998i . We note, however, that in spite of the general success 
of the Shakura & Sunyaev a-parametrization of turbulence, in 
particular for hot disks around compact objects, it is known that 
the this model does not fit very well all of the observed disks 
of TTauri stars (i.e. cool disks), in which a T(r) ~ r" 1 ' 2 tem- 
perature profile indicated (Beckwith et al. H990i . 

Finally, we supply the expression for the disk-half thickness 
from Shakura & Sunyaev for the case of a circum-planetary 



3/20 



h(r) = 6.5 x 10 6 cma 



-1/10 




.(A.5) 



Compared to the circum-stellar accretion disk, the circum- 
planetary is expected to somewhat thinner. 
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